arXiv:1505.00938vl [cond-mat.mes-hall] 5 May 2015 


Surface excitations in the modelling of electron transport for electron- 
beam-induced deposition experiments 

Francesc Salvat-Pujol* 1 , Roser Valenti 1 and Wolfgang S. Werner 2 

Address: 1 Institut fur Theoretische Physik, Goethe-Universitat Frankfurt, Max-von-Laue-StraBe 1, 
60438 Frankfurt am Main, Germany and 2 Institut fur Angewandte Physik, Technische Universitat 
Wien, Wiedner HauptstraBe 8-10/134, 1040 Wien, Austria 

Email: Francesc Salvat-Pujol - salvat-pujol@itp.uni-frankfurt.de 
* Corresponding author 

Abstract 

The aim of the present overview article is to raise awareness of an essential aspect that is usually 
not accounted for in the modelling of electron transport for focused-electron-beam-induced depo¬ 
sition (FEBID) of nanostructures: surface excitations are on the one hand responsible for a size¬ 
able fraction of the intensity in reflection-electron-energy-loss spectra for primary electron ener¬ 
gies of up to a few keV and, on the other hand, they play a key role in the emission of secondary 
electrons from solids, regardless of the primary energy. In this overview work we present a gen¬ 
eral perspective of recent works on the subject of surface excitations and on low-energy electron 
transport, highlighting the most relevant aspects for the modelling of electron transport in FEBID 
simulations. 
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Introduction 


An accurate modelling of the energy losses of electrons traversing a solid surface is instrumental 
for a quantitative understanding of a series of techniques exploiting transmitted, reflected, or emit¬ 
ted electrons, including a number of spectroscopies (electron-energy-loss spectroscopy, x-ray pho¬ 
toelectron spectroscopy, and Auger-electron spectroscopy), electron microscopy, and the focused- 
electron-beam-induced deposition (FEBID) of nanostructures, on which we focus here. This tech¬ 
nique employs beams of focussed keV-electrons to trigger and steer the growth of nanostructures 
with tunable electronic and magnetic properties from molecules of organometallic precursor gases 
[|T|| adsorbed on a substrate J2). It has been shown that, for irradiation with electrons of 1-5 keV, 
both the incoming primary electrons and the emitted secondary electrons influence the growth of 
the nanostructure, the latter electrons being responsible for the lateral resolution [[3]|. 

In the modelling of electron transport for FEBID [[2)13, electron stopping is described on the ba¬ 
sis of properties that are applicable in the bulk of the material. However, electrons traversing a 
solid interface additionally excite surface modes, an energy-loss channel that amounts to a size¬ 
able fraction of the energy-loss spectrum for electrons of up to a few keV. The existence of surface 
excitations was predicted in the late 1950s JS); experimental evidence was obtained shortly there¬ 
after 11911 1011 . In order to model surface excitations in electron spectroscopies, several models have 
been developed to date HI 1113 111 . Various approaches are considered, often with underlying simpli¬ 
fications, evoked on physical or technical grounds, in the interest of making calculations feasible 
in a finite time. In order to derive a distribution of energy losses of charged projectiles moving in 
the vicinity of the surface, some of the models cited above rely on the semiclassical dielectric for¬ 
malism, whereas others adopt a many-body formalism. Both approaches have been shown to yield 
results in equivalently good agreement ll32ll with experimental data. 

In what follows we briefly review the stopping of charged projectiles in the vicinity of a solid sur¬ 
face, along the lines of lf3TI . which will be referenced for further details. We summarize a series 
of rules which characterize the behavior of the probability for surface excitations and we briefly 
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review a practical model for the emission of secondary electrons. Relevant aspects to FEBID mod¬ 
elling will be highlighted. 

Inelastic collisions in the bulk of the material 

Energy losses of a charged projectile moving in a solid can be described accurately within the 
semiclassical dielectric formalism. In this approach, one assumes that the presence of the charged 
projectile disturbs the equilibrium charge density of the solid, which becomes polarized and, thus, 
an electric field is induced at all points of space. The force acting on the charged projectile due to 
the induced electric field is assumed to be the agent responsible for its (electronic) stopping. In or¬ 
der to derive physical quantities that describe the stopping, it is now a matter of calculating first 
the induced electric field and, from it, the so-called stopping power, defined as the variation of the 
kinetic energy of the projectile per unit path length. Once an expression for the stopping power is 
derived, one can identify from it an expression for the distribution of energy losses per unit path 
length, the basic quantity that is needed to describe energy losses in a detailed Monte-Carlo simula¬ 
tion of electron transport. In this section we briefly outline the basic steps of these calculations and 
highlight the underlying assumptions. Further details can be found in the cited works. 

The starting point of the calculation is the dielectric function e(q, co ) of the material, where q and 
co are the respective Fourier-conjugate variables of the position, r, and the time, t. In practice one 
typically has data available for e(a)), be it from optical data obtained experimentally f[33l or from 
theoretical calculations, e. g., via density-functional theory calculations H33H35H . An co-dependent 
dielectric function is sufficient to describe the response of the medium to a spatially homogeneous 
perturbation, such as that of an incoming photon. However, for incoming charged projectiles the 
perturbation is strongly dependent on the spatial coordinates, so that a ^-dependent dielectric func¬ 
tions is required. Physically reliable models are built on the basis of the (q, co )-dependent dielec¬ 
tric function for the homogeneous electron gas H36H381 or on the basis of a simple superposition of 
Drude-Findhard oscillators If33ll . 

Assuming a projectile that moves with a velocity v along a trajectory r = \t, one can conveniently 
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solve the Maxwell equations in Fourier space to obtain the following expression for the induced 
electric field OH 


E in «i(q,ffl) = —i47T-^-p(q, a>) 
<T 


£(q,w) 


-1 


( 1 ) 


where p(q, (o) is the Fourier transform of the projectile charge density p(r,f) = Zo<?<5(r — \t), 
where Zq is the projectile charge in units of the modulus of the electron charge, e, and v is the ve¬ 
locity of the projectile. To obtain this expression, the following approximations were considered: 

(1) the Coulomb gauge was adopted and the contribution from the vector potential was neglected 

(2) the dielectric displacement field was assumed to be proportional to the electric field in Fourier 
space (linear response). The first approximation restricts the validity of the calculation to non- 
relativistic projectiles (the calculation with the full electric field for relativistic projectiles is also 
feasible ll39l ). whereas the second approximation can be seen to be formally equivalent to a first- 
order Born approximation in perturbation theory, imposing a lower bound to the domain of validity 
of the calculation H22I401 . which for practical purposes is above 100 eV. 

The stopping power S is obtained as the variation of the kinetic energy per unit path length, 
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d<§ 

di 


y-E ind (iy)| r=vf , 


( 2 ) 


where $ is the kinetic energy of the projectile and s — vt is the path length. Combining Eqs. (1) 
and (2) one obtains 
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Up to this point the stopping of the projectile is treated as a continuous phenomenon, whereas in 
reality charged projectiles lose energy and are deflected in the course of individual inelastic colli¬ 
sions. The so-called semiclassical approximation consists in assigning to ft q and ftO) the meaning 
of a momentum transfer from the projectile to the medium and of an energy loss of the projectile, 


4 











respectively. Atomic units (h = m e — e — 1) will be used below. Now that these variables have a 
well-defined physical meaning, the corresponding integrals must be restricted to the kinematically 
allowed domain. 
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where 


q± = \/l.S ± \jl{$ — to) 


(5) 


are the minimum (-) and maximum (+) allowed momentum transfers allowed by the energy and 
momentum conservation laws. 

Equation (4) can be understood as the average energy loss per unit path length dictated by a distri¬ 
bution of energy losses per unit path length, d/i/dft): 
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The quantity d/i/dco is known as the differential inelastic inverse mean free path (DIIMFP), explic¬ 
itly given by 
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Note that the DIIMFP is a function of the energy loss for the given velocity of the projectile. The 
integral of the DIIMFP over all allowed energy losses gives the inelastic inverse mean free path 
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The latter two quantities are the necessary quantities for a detailed Monte Carlo simulation of elec- 
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tron transport (see section “Monte-Carlo simulation of electron energy-loss spectra”), a method 
that has been successfully used in the last decades. 


Inelastic collisions in the vicinity of a planar surface 

The scheme outlined in the previous section to describe inelastic interactions of charged projec¬ 
tiles in solids gives a good account of inelastic collisions in the bulk of the solid. However, projec¬ 
tiles impinging and emerging from a solid additionally cross a planar interface to vacuum (or an¬ 
other solid) that is not explicitly accounted for. The existence of a plane surface imposes additional 
boundary conditions on the electric field H411421 . 

Several approaches exist in the literature to solve the Maxwell equations with these boundary con¬ 
ditions for the stopping problem: some consider the dielectric function of a semi-infinite medium 
@31, and others (preferred in the electron-spectroscopy community) rely on a method which allows 
one to work with bulk dielectric functions, the method of image charges, also known as the method 
of extended pseudomedia. The method consists in rephrasing the semi-infinite-geometry problem 
as the sum of two infinite-geometry problems, supplied with a series of fictitious charges that are 
determined in terms of known quantities by imposing the boundary conditions at the interface. 

The resulting induced electric field has a more complex expression than in the bulk case. Never¬ 
theless, it can be expressed as the sum of one contribution arising from a charge density induced in 
the bulk of the material and another one arising from a charge density induced at the surface of the 
material. 

The DIIMFP resulting from the induced electric field becomes more complicated, with two addi¬ 
tional parametric dependencies: (1) on the depth coordinate with respect to the surface and (2) on 
the surface crossing angle with respect to the surface normal. Several models exist with varying ap¬ 
proximations H23I24I30I31L the effect of which was scrutinized iPTl . Regardless of the details of 
the models, they all yield a number of consistent general features and trends of the surface excita¬ 
tion probability: 

• Surface energy losses can be undergone by the charged projectile on either side of the inter- 
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face, at the solid side or at the vacuum side (!). Indeed, a surface charge can be induced re¬ 
gardless of the side at which the projectile is moving on and, thus, a charged projectile mov¬ 
ing on the vacuum side of the interface can also undergo energy losses. It has been recently 
shown that, in reflection-electron-energy-loss spectra, surface losses on the vacuum side of 
the interface account for a large fraction of the surface-excitation intensity, often more than 
half of it EM- 

• The probability for an electron that crosses a surface to undergo a surface excitation is, to a 
first approximation ll45ll . proportional to the surface dwell time l rs_/ l/(V^cos0), where S 
is the projectile energy and 0 is the surface crossing angle with respect to the surface normal. 
The energy dependency implies that, in practice, surface excitations are relevant for elec¬ 
tron energies up to a few keV. Additional structure to the aforementioned angular behavior 

is predicted for scattering geometries coinciding with deep minima of the differential elastic 
scattering cross section: minor deflections in the course of an inelastic collision lead to an 
effective scattering geometry with enhanced elastic scattering and therefore higher detection 
probability lf44il . 

• The DIIMFP for energy losses of charged projectiles impinging on a surface differs from 
the DIIMFP for the conjugate emerging direction. This effect, known as in-out asymmetry 
in surface energy-losses, has been long predicted but only recently observed experimentally 
[f46ll . In-out differences are most accentuated for surface-crossing directions close to the sur¬ 
face normal and for high kinetic energies (~ 1 keV). 

Monte-Carlo simulation of electron energy-loss spectra 

The electron-transport problem in a solid is described in terms of a Boltzmann-type transport equa¬ 
tion. A practical method for solving the problem is provided by Monte-Carlo simulation, which 
consists in sampling an ensemble of trajectories undergoing collisions of the relevant types as dic- 
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Vacuum 



~±15 A 


Figure 1: Example incoming trajectories (dotted lines) in the surface-scattering zone (typically 
±15 A about the surface), undergoing (a) a surface energy loss in the medium side, (b) a surface 
energy loss in the vacuum side, (c) a bulk energy loss, indicated by the filled circles. 


tated by a given set of interaction cross sections. A statistical average of the desired observable is 
performed over the sampled trajectories to the selected precision [|47l . 

In the energy range between 100 eV and a few keV the relevant interaction mechanisms of elec¬ 
trons with the solid are elastic collisions with the atoms and inelastic collisions with typically 
weakly bound electrons in the solid. Elastic scattering can be accurately described by means of a 
differential cross section for elastic scattering (DCES), which can be systematically calculated by 
means of partial-wave calculations H48I49I . Inelastic scattering is accounted for by the DIIMFPs 
described above. Monte-Carlo simulations of electron transport (bulk losses only) for typical ge¬ 
ometries in FEBID experiments have been previously considered Il2l71 . The inclusion of surface 
excitations implies a modification of the sampling algorithm in the vicinity of the surface (typi¬ 
cally ±15A about the surface), as schematically shown in Fig. [T] Technical details on the imple¬ 
mentation of the algorithm for the simulation of surface energy losses can be found elsewhere 
in great detail 001311501 . Here the focus is on the effect of surface excitations on the reflection- 
electron-energy-loss spectrum (REELS). To this effect, Fig. [2] compares the REELS of Si (left) and 
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Cu (right) under bombardment with 1-keV electrons impinging perpendicularly on the sample; all 
reflected electrons are collected. The simulation geometry is depicted in Fig. [3] The materials are 
chosen as representative substrate (Si) and deposit (Cu) materials. The solid red curves (dashed 
blue curves) in Fig. [2] correspond to REELS simulated without (with) the inclusion of surface exci¬ 
tations. We observe that even for a primary energy of 1 keV surface excitations account for (1) ad¬ 
ditional features, i. e. the excitation of surface plasmons, in the low-energy-loss part of the REELS 
that are not accounted for by a bulk-only description of the energy losses of charged projectiles in 
the material and (2) a sizeable fraction of the intensity in the first few tens of eV of energy losses, 
about 20% of the intensity in the case of Si and 15% of the intensity in the case of Cu. Although 
the relative importance of surface excitations is enhanced for lower energies, their effect is no¬ 
ticeable even in the 1-keV domain. Thus, the inclusion of surface excitations in the modelling of 
electron-transport is expected to give a yet more quantitative description of FEBID processes at 
and below the 1-keV primary-energy domain. 

Si, £7 = 1 keV Cu, E — 1 keV 




Figure 2: Comparison of reflection-electron-energy-loss spectra (REELS) of Si (left) and Cu 
(right) under bombardment with 1-keV electrons at normal incidence, without (red solid curves) 
and with (blue dashed curves) an account of surface excitations. All backscattered electrons are 
collected. 
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Figure 3: Simulation geometry: 1-keV electrons impinge normally onto the material (Si or Cu); all 
backscattered electrons are collected. 


Secondary-electron emission 

Energy losses of the charged projectile can lead to the ejection of loosely bound electrons of the 
solid, which emerge as secondary electrons (SE). The majority of these SE are of relatively low en¬ 
ergies (< 50 eV). These energies are well below the domain of validity of the elastic and inelastic 
interaction cross sections available in the literature, which has been a limitation for progress in the 
field. Electron coincidence measurements H51H54H have supplied a wealth of valuable information. 
Recently, coincidence measurements of correlated electron pairs emitted from solids (Al, Si, Ag) 
under electron bombardment have been measured, providing a double-differential SE yield, differ¬ 
ential with respect to the energy loss of the primary electron and with respect to the energy (or the 
time of flight) of the emitted secondary electron 11551 . These experimental data are displayed for Si 
under 100-eV electron bombardment in the lower panel of Fig. |4]as a bird’s-eye-view plot (only the 
shape and relative intensities of the spectrum are of relevance here, hence the missing units in the 
linear color scale, where black is the null point and white is the maximum attained value). The hor¬ 
izontal white lines indicate the corresponding times of flight for electrons with 0 eV (accelerating 
grids were used), 50 eV, and 100 eV. See Il55l for the experimental details. The plot can be read as 
the (time-of-flight) spectrum of secondary electrons emitted as a result of different energy losses of 
the impinging electron (to be read at the abscissae). The upper panel of Fig. [4] displays the REELS 
of 100 eV from Si, where the energy-loss peaks corresponding to the excitation of one surface plas- 
mon, one bulk plasmon, and two surface plasmons are indicated by vertical red dashed lines and 
labeled, respectively, Is, lb, 2s as a guideline for the abscissae scale in the other plots of the figure. 
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Figure 4: (Upper panel) Reflection-electron-energy-loss spectrum (REELS) of Si under 100 eV 
bombardment (see If55l for experimental details). (Lower panel) (e,2e)-coincidence spectrum of 
secondary electrons emitted in coincidence with energy losses (SE2ELCS) of 100-eV electrons 
backscattered from Si. (Middle panel) Monte-Carlo simulation of the SE2ELCS measurement 
without accounting for surface energy losses. The vertical dashed lines in red indicate energy 
losses corresponding to the excitation of one surface, one bulk, and two surface plasmons. The 
horizontal solid white lines indicate the times of flight corresponding to electrons with 0 eV (ac¬ 
celerating grids were used), 10 eV, and 50 eV. 


The coincidence data (e. g ., lower panel of Fig. [4]) provide on the one hand very detailed insight 
into the mechanisms responsible for SE emission and, on the other hand, provide a benchmark 
against which models for SE emission and low-energy electron transport in general can be tested. 
The transport models, in turn, aid in the interpretation of the data, as discussed below. The Monte- 
Carlo simulation briefly outlined above was extended to include the generation and the transport 
of the secondary electrons and to simulate the electron-coincidence measurement on the basis of 
a simple model for SE emission: every time that the primary electron undergoes an energy loss, 
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Figure 5: Same as Fig. [4] with the inclusion of surface excitations in the modelling of electron 
transport through the solid-vacuum interface. 

a SE trajectory is started with the energy loss as an initial energy (see If55ll for the simulation de¬ 
tails). Having the experimental data as a guideline, the interaction cross sections described above 
were used down to 1 eY (knowing that this is well below the domain where they are formally ap¬ 
plicable) as a first approximation. Simulations were first carried out using bulk energy-loss DI- 
IMFPs exclusively. The resulting spectrum is displayed in the middle panel of Fig. |4j It is clear 
that these simulated peaks do not reproduce the onset of the experimentally observed peaks. Only 
after the inclusion of surface excitations, both for the incoming primary electrons, for the backscat- 
tered electrons, and for the emitted secondary electrons, is good agreement between simulations 
and measurements found, as shown in Fig. [5] The Monte Carlo simulations further allow one to 
discern the processes that give rise to the different regions of the coincidence spectrum lf55ll . 

Thus, it was found that any realistic model of SE emission and low-energy electron transport near 
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solid surfaces must account for surface excitations. This conclusion has strong implications on 
the emission depth from which SE are emitted: if secondary electrons undergo additional energy 
losses on their way out of the solid, the average SE-emission depth becomes much shallower than 
one would assume on the basis of a model based only on bulk properties. The predicted number 
of emitted SE can also differ appreciably with respect to a bulk-only simulation. Furthermore, the 
energies of the SE are also modified by the presence of additional surface energy-loss channels. 

Conclusions 

In light of the presented richness in the behavior of surface excitations and their effect on both the 
energy losses of the impinging electrons and on the emission of secondary electrons, it is to be ex¬ 
pected that their inclusion in the modelling of electron transport for FEBID will yield a more de¬ 
tailed description of the role played by both the primary electrons and the emitted secondary elec¬ 
trons in the growth process. It should be noted that, while surface excitations are relevant for pri¬ 
mary electrons with energies up to 1-2 keV, they are essential ingredients for the modelling of slow 
secondary electrons regardless of the energy of the primary electron responsible for their emission. 
The previous considerations suggest that the inclusion of surface excitations in the electron- 
transport model employed to investigate FEBID experiments might lead to noticeable effects. On 
the one hand, more primary electrons are backscattered compared to the case without surface exci¬ 
tations (see Fig. 2), so that an increase in the simulated deposition rate might be expected (at least 
for primary energies in the 1-2 keV regime and below). On the other hand, more slow (<= 50 eV) 
secondary electrons will be available from the decay of surface plasmons lf56ll excited by either the 
incoming electrons or the outgoing electrons (backscattered electrons or emitted secondary elec¬ 
trons). This should also contribute to an increase of the simulated deposition rate and additionally 
lead to an enhancement of the FEBID proximity effect. 
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